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Abstract 



It is proved that the cosmological density perturbation is associated with 
a pecuhar velocity field. This allows a simple formulation of cosmological 
perturbation theory, which works entirely with quasi-Newtonian fluid flow 
equations. As an illustration, the large scale cosmic microwave background 
anisotropy (Sachs- Wolfe effect) is calculated without any reference to the met- 
ric perturbation. In addition, assuming the usual adiabatic initial condition 
on the density perturbation, we show that the dipole of the anisotropy mea- 
sures our peculiar velocity relative to the average peculiar velocity within the 
last scattering surface of the microwave background, thus defining its frame 
of reference. 
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The interplay between particle physics and cosmology, generally known as particle cos- 
mology, is particularly fruitful at the present time, and is likely to play a decisive role in 
both fields during the next few years. For some time, the main focus has has been on the 
origin of structure in the universe. In this context, one is seeking a huge body of data on 
galaxies and galaxy clusters, as well as the cosmic microwave background (cmb) anisotropy. 
The first observation [|I| in 1992 of the intrinsic part of this isotropy has been one reason 
for for a dramatic surge of interest in structure formation, and the ever improving quality 
of data on galaxies and clusters has been another. 

The simplest hypothesis concerning the origin of large scale structure is to invoke a 
primeval density perturbation (plus perhaps a gravitational wave contribution to the cmb 
anisotropy), which originates during inflation as a vacuum fluctuation. If this hypothesis 
is correct, the observations mentioned above provide a window on the nature of the funda- 
mental interactions at an energy scale which is many orders of magnitude bigger than those 
which can be explored at colliders. 

In order to explore the hypothesis, one needs to formulate a relativistic theory of cosmo- 
logical perturbations. At present three rival formulations exist in the literature. The 'metric 
perturbation' formulation, originated by Lifshitz in 1946 0, considers the components of 
the metric tensor, which are related to the components of the energy-momentum tensor in 
either the perfect fluid approximation, or including the effects of particle diffusion and free- 
streaming 1^. The equations in this approach are complicated, though they are perfectly 
serviceable and have in fact been used in most of the decisive comparisons of theory with 
observation. 

Second, there is the 'gauge invariant' approach initiated by Bardeen in 1980 ^ and 
extended by various authors to include diffusion and free- streaming |^. It starts with the 
metric perturbation approach, but allows one to derive a much simpler set of equations for 
matter perturbation variables which do not explicitly involve the metric perturbation. The 
only problem is that the physical interpretation of the variables used is somewhat obscure, 
as is that of the mathematical manipulations through which the equations are derived. 

Finally there is an approach pioneered by Hawking in 1966 but developed mostly 
in the last decade [[^-[T2[| which is based on the 'covariant' fluid flow approach to general 
relativity [|1^]. Here, one bypasses the metric tensor to work directly with the curvature 
tensor, which allows a simple treatment of perturbations, even permitting a coordinate- 
free, gauge invariant description ||,|I^,|rT[] . This approach has been developed for a single 
perfect fluid [^-|T0| and for several uncoupled perfect fluids [0,|ll[, and in both cases it has 



been shown to provide a simple, physically transparent, derivation of the main results of 
the 'gauge invariant' approach. Although all of the 'gauge invariant' variables have been 



recently identifled with physical quantities [10,11|, no relation has yet been provided with 



one of the most important observables, i.e. with the peculiar velocity fleld; in addition, 
diffusion and free streaming remain to be treated. 

Something which has hitherto generally being lacking in relativistic perturbation theory 
is a straightforward procedure for going to the Newtonian limit. Moreover, the concept of 
peculiar velocity, which plays a vital role in the Newtonian case, has not so far even been 
mentioned in the relativistic context. 

In this letter we show, using the 'covariant' approach, how to extend the peculiar velocity 
(PV) concept to the relativistic regime. As an application of the concept, we clarify the 
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relation between the PV field and the cosmic microwave background (CMB) anisotropy, 
deriving the Sachs Wolfe effect in a straightforward way, without any mention to metric 
perturbations. Furthermore, we show that the dipole of the anisotropy measures our PV 
with respect to the average PV within the last scattering surface of the CMB, thus defining 
for the first time a cosmic frame of reference. We end by relating our peculiar velocity to 
one of the 'gauge invariant' variables, and pointing to the wider implications of our work. 

This letter is necessarily brief with emphasis on the Newtonian analogy: full results, 
as well as a clarification of the link between the treatment here and the other possible 
approaches mentioned before, will be reported elsewhere 

The universe is regarded as a fluid, taken here to be perfect, and a comoving observer 
measures by definition zero momentum density, energy density p and pressure p. At a given 
spacetime point one can use locally inertial coordinates in which the fiuid is instanta- 
neously at rest. Then |T3 the fiuid 4- velocity at (infinitesimally) nearby spacetime points 



defines the acceleration Oj = dui/dy^ of the fiuid, and its velocity gradient Uij = duj/dy"^. 
The velocity gradient can be uniquely decomposed into an antisymmetric vorticity Uij, a 
symmetric traceless shear ajj, and a locally defined Hubble parameter if, 

Uij = H6ij + cxij + ujij . (1) 

The unperturbed universe is isotropic around each comoving observer, in particular the 
acceleration, shear and vorticity all vanish. Since the vorticity vanishes there exist hyper- 
surfaces of simultaneity, orthogonal to the fiuid worldlines. Isotropy about every comoving 
observer implies that these hypersurfaces are homogeneous; on them p, p and H are position- 
independent, and so is the proper time t along a fiuid worldline, starting from one of the 
hypersurfaces. Making the assumption of critical density = (87rG/3)p, the hypersurfaces 
are fiat and one can define on them coordinates Xj which are fixed along each fiuid worldline 
(comoving), such that Cartesian coordinates (the generalisation to non-critical 

density is straightforward 0]). The scale factor, defined by aH = da/dt, is normalised to 
1 at the present epoch. 

Newtonian physics is valid after matter domination, on scales much smaller than the 
Hubble distance H^^. In this case there is a well defined fluid velocity u, and a PV field v 
related to it by 

u(r) - u(0) = /fr + v(r) - v(0) , (2) 

where r is the displacement from our position (or that of any other comoving observer), and 
H the average expansion. An equivalent statement in terms of the velocity gradient is 

Uij = 6ijH + diVj . (3) 

These expressions define v up to a constant, which can be chosen so that the average of v 
vanishes. Making the decomposition |16| v = v^ + v""", where the transverse part v""" satisfies 



divJ = and the longitudinal part is of the form = Vt/'^,, the comparison with (|T]) gives 



1 



^ij = ^{div- - djvl) , (4) 
5H = ^V^^. , (5) 
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where 6H is the perturbation in the expansion, while the shear aij does not carry any new 
information. Defining the pecuhar gravitational potential ip by 

A7iG6p = V'^ip , (6) 

the Newtonian hydrodynamical equations show that it is constant (dropping the decaying 
mode), and that 

vL(x,t) = -tV^(x). (7) 

Now consider the relativistic regime. In the presence of perturbations, the special coor- 
dinates t and Xi do not exist, but any first order approximation to them may be used to 
describe the perturbations, with an error of only second order. First order perturbations 
'live' in the unperturbed universe. We use the notation di = d/dri = a~^d/dxi, and let an 
overdot denote differentiation with respect to t at fixed x. Writing each perturbation /(x, t) 
as a Fourier series in a comoving box much bigger than the region of interest, modes with 
different wave- vector k/a decouple in the equations. A mode is said to 'enter the horizon' 
when aH/k falls below 1, and for scales > 200 Mpc this occurs after matter-domination. 

Because they are time dependent, perturbations in p, p and H have to be defined with 
respect to a slicing of spacetime into spacelike hypersurfaces which become homoge- 
neous in the unperturbed limit. On each hypersurface we split p into an average plus a 
perturbation, p = p + 6p, and similarly for p, H. 

We are now going to show that the slicing can be chosen so that there exists a relativistic 
generalization of Eq. (0) which defines a unique longitudinal PV = V?/'^. This relativistic 
generalization of the Newtonian concept of PV is the central result of this letter. 

If the vorticity does not vanish there are no hypersurfaces orthogonal to the fiuid world- 
lines, but Eq. §) defines a unique v . Then we can define modified fiuid worldlines 



which at each spacetime point have velocity — v"'" with respect to a comoving observer. 
These worldlines have zero vorticity, so there exist hypersurfaces orthogonal to them, called 
comoving hypersurfaces |T^, on which we define Sp, 6p and 6H. 

Equations involving only ujij, aij, 6p, 6p and 6H can be derived from the Einstein 



field equation, together with the Ricci identity acting on the 4- velocity u^. One of them is 
the constraint equation 

2di6H + HvJ = djicxij + LUij) , (8) 

(The left hand side is the derivative perpendicular to the fiuid worldlines, whereas the 
first term alone is the derivative perpendicular to the modified worldlines, i.e., within the 
comoving hypersurfaces). As we now show, this equation is equivalent to Eq. (^ in the 
Newtonian limit, and gives the desired generalisation of it in the relativistic regime. We 
start with the fact that any traceless symmetric tensor field may be decomposed uniquely 
p!6| in the form 

(^ij = ^idiwj + djwj) + VijXv + ol , (9) 
where Vjj = didj — ^^jj-V^, and w'^ and aj^ are transverse, diwj = and diajj = 0. Also, 



any scalar field 6H may |16| be written uniquely in the form Eq. (^. Then, using Eqs. 
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(I) and @, the longitudinal part of Eq. (||) is Xv = "ipv, and its transverse part is in Fourier 
space 



w 



T 
V . 



(10) 



Calculating Suij = 6HSij + aij + Uij one therefore finds 



6u 



-{divj - djvj) + -{diwj + djwj) . 



(11) 



Thus, in the relativistic regime there is a well defined longitudinal PV v^, but two 
different 'transverse peculiar velocities' v""" and w""", related by Eq. (0). The difference is a 
purely relativistic effect (the dragging of inertial frames) and one verifies from Eq. (p!OD that 
it disappears in the Newtonian regime aH/k <^ 1. 

The symmetric, traceless, transverse term ajj directly represents the effect on matter 
motion of gravitational waves, a degree of freedom absent in Newtonian physics. 

The other equations involving only Uij, aij, Sp, Sp and 6H are propagation equations. 



that in the Newtonian limit are equivalent to the usual fluid flow equations. One of them ||T0[ 
involves only the vorticity Uij, and shows that it decays like {p + p)~ 
means that it is negligible. The others are 0, 



which presumably 



{SpJ= -3{p + p)SH -3HSp 
(SH) = -2H6H 



4:7iG , 1 \/^5p 



3 p + p 



(12) 
(13) 



After matter domination they lead to the Newtonian expressions Eqs. (^) and (|^), even on 
scales aH/ ^ 1 where Newtonian physics does not apply. 

The constraint equation and the three evolution equations completely determine the 
evolution of the velocity gradient, except for the gravitational wave contribution ajj. This is 



related to the traceless transverse part hj^ of the spatial metric perturbation by af^ 



T0| , which in turn satisfies pO 



1/,TT 



TT 



2 7, TT 



(14) 



Well before horizon entry, each Fourier component hjj{t, k) has some constant value Aij(k) 
(ignoring a decaying mode). Well afterwards, it oscillates as a standing wave with am- 
plitude decreasing like a~^, and wavenumber and angular frequency k/a. For the scales 
> lOOMpc which enter the horizon after matter domination (when a^^ oc {aHY) its 
contribution to the velocity gradient is therefore 5Mjj(grav) ~ AijH^a/k. 

Let us ask how significant this contribution is compared with the density perturbation 



,-1 



contribution Suij{dens) 



From Eq. (0), 

6uij (giav) 

r 

6uij{dens) 



(15) 
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In general there is no reason why the ratio should not be large at horizon entry, in which 
case the gravitational waves could still be important well after horizon entry. However, if 



ip and Aij both originate as a vacuum fluctuation during inflation, then [13,20,5^ Ay <^ 



(on the scales k > 100 Mpc that we are considering.) This implies that on the scale 
~ 100 Mpc, 6uij{grav) is < 10~^6uij{deias) at the present epoch, amply justifying its 
neglect when is deduced from large scale galaxy surveys. 

To demonstrate the utility of Eq. (Ill]), we now use it ||12| to calculate the Sachs- Wolfe 



effect |22|, which is normally written in terms of the metric perturbation, and describes 
part of the anisotropy of the CMB. This anisotropy is defined as the angular variation in 
the intensity of the background at fixed wavelength A, and is usually specified by giving 
the equivalent variation A(e) = AT/T in the temperature T of the blackbody distribution, 
where e is a unit vector pointing in the direction of observation. Some anisotropy Aem(e) is 
already present on a comoving hypersurface just after last scattering, and the Sachs- Wolfe 
effect describes the additional anisotropy acquired on the journey towards us, to first order 
in the perturbations. The Sachs- Wolfe effect is due entirely to the anisotropy of the redshift 
of the CMB, as the photons path can be taken as unperturbed in first order. 

Consider a photon passing a succession of comoving observers. Its trajectory is dr/dt = 
— e and between nearby observers its Doppler shift is 

d\ ^ da ^ ^ 

— = CiCjUijdr = — — + eiCjOUijdr , (lb) 

A a 

where the first term is due to the average expansion, and the second is due to the relative 
PV of the observers. Integrating this expression gives the redshift of radiation received by 
us, which was emitted from a distant comoving source. The unperturbed result is A/Agm = 
1/aem, and the first order perturbation gives the Sachs- Wolfe effect 

A(e) - Acm(e) = / eiej5uij{yi,t)a{t)dx , (17) 

<J 

where Xgm — '^Hq^ is the coordinate distance of the last scattering surface, and the integra- 
tion is along the photon trajectory 

^5 = 3 f^-iV (18) 



The symmetric part of Eq. (11) gives the various contributions to the Sachs- Wolfe ef- 
fect. The transverse PV is presumably absent, and we have nothing new to say about the 
gravitational wave contribution. This leaves v'", which after integrating by parts and using 
Eqs. ([T|) and (0) gives 

A(e) = Ae^(e) + ^[z^(xe^) - ^/.(O)] + 

e- [v(0,to) - v(Xen„tem)] . (19) 

Here Xem = a;eme, and for clarity v'" is denoted by v. 

A better expression follows if one uses the divergence theorem and Eq. (|l^) to project 
out the dipole part of '?/'(xem)/3. One finds that it is equal to (v(tejn) — v(to)) where () 
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denotes the average within the last scattering surface x = Xem- Defining v = v — (v) this 
gives 

A(e) = Acm(e) + ^ [V^(Xem)]i>i + 

e ■ [v(0, to) - v(Xem, tem)] • (20) 

On angular scales > 1°, corresponding to linear scales at last scattering which are outside 
the horizon, v(xein,tcm) is negligible compared with ip, and with the usual adiabatic initial 
condition so is Aem(e) ~ |5p(xcm5 4m)/p(^em) Dropping the monopole one then 

obtains 

A(e) = ^[V^(xem)],>i + e-vo. (21) 

The last term in this expression is the dipole, and defines for the first time the rest frame 
of the CMB (in linear theory with the adiabatic initial condition). This is simply the frame 
at rest with respect to the average peculiar velocity of everything within the last scattering 
surface. Of course this result corresponds roughly to the intuitive idea that a distant density 
perturbation will affect the matter and the CMB almost equally ||23|| . 

Although the above application of the peculiar velocity concept is amusing, it is far from 
being the most important aspect of our work and we end by briefly summarizing the wider 
picture, dropping for simplicity the rotational and gravitational wave modes. 



Here we have been using the 'covariant' |]T2| fluid flow approach to cosmological pertur- 



bation theory as opposed to the 'metric' [Q] and 'gauge invariant' [|| approaches, showing 
how close is our treatment to the Newtonian one. The link between the two latter has 
been deeply analyzed 0], while the relations between the main variables of the 'covariant' 
approach and those of the 'gauge invariant' treatment were given [0,|ll|], but not yet fully 
exploited to treat diffusion and free-streaming. 

Our new point of contact with the 'gauge invariant' approach is the fact that the Bardeen 
1^ variable v^^^ is related to our peculiar velocity potential ipv- In fact, Bardeen pointed 
out in that original reference (Eqs. (3.11) plus the one before Eq. (3.12)) that the shear of 
the velocity fleld associated with v^^^ is just the shear of the comoving worldlines. What 
we have shown here is that, in a comoving slicing |T^ and a critical density (flat) universe. 



there is a unique quantity = 'Vipy, which gives both the shear of the comoving worldlines 
and the perturbation in their expansion through the relation 36H = V.v. It follows that, 
under these conditions, the velocity field associated with vf'^ is precisely the pecuhar velocity 
defined in the present paper. 



Adding this connection to the already known one |TO|,|Tl|, it should now be possible to 
write and simply derive all of the equations of the gauge invariant approach following a 
covariant treatment. The result will be a complete description of cosmological perturbation 
theory, simpler than either of the alternatives and treating the Newtonian and relativistic 
regimes on essentially the same footing. 

While this manuscript was in preparation, we became aware that a related derivation 
of the Sachs- Wolfe effect was under study [Q. We are grateful to P. Coles, F. Lucchin, R. 



Scaramella for useful comments. MB thanks SERC (UK) and MURST (Italy) for financial 
support. 
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